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PREFACE 


The  work  reported  herein  was  performed  in  the  Hydraulics  Laboratory  (HL) 
of  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  as  a  part  of  the 
overall  investigation  to  predict  the  evolution  of  the  Atchafalaya  Bay  delta. 
The  study  design  (Phase  I  of  the  study)  was  authorized  by  the  US  Army  Engineer 
District,  New  Orleans  (LMN),  on  18  July  1977.  The  implementation  of  the  study 
plan  (Phase  II)  was  authorized  by  LMN  on  21  May  1979.  This  report  presents 
work  done  under  Phase  II,  in  support  of  numerical  modeling  of  the  delta 
evolution. 

This  study  was  conducted  under  the  direction  of  Messrs.  H.  B.  Simmons 
and  F.  A.  Herrmann,  Jr.,  former  and  present  Chiefs,  HL;  R.  A.  Sager,  Assistant 
Chief,  HL;  W.  H.  McAnally,  Chief,  Estuaries  Division;  G.  M.  Fisackerly,  Chief, 
Estuarine  Processes  Branch;  and  R.  A.  Boland,  former  Chief,  Hydrodynamics 
Branch.  The  plan  of  study  of  which  this  task  is  one  part  was  developed  by 
Messrs.  McAnally  and  Samuel  B.  Heltzel,  Estuarine  Engineering  Branch.  This 
study  was  performed  by  Mr.  A.  M.  Teeter,  Estuarine  Processes  Branch,  assisted 
by  Messrs.  B„  G.  Moore  and  D.  M.  White,  Estuaries  Division.  Gratitude  is 
extended  to  Dr.  R.  Ariathurai  and  Mr.  McAnally  for  consultation.  Consultants 
to  the  project  were  Mr.  L,  R.  Beard,  Dr.  C.  R.  Kolb,  Dr.  R.  B.  Krone,  and 
Mr.  F.  B.  Toffaletti.  This  report  was  written  by  Messrs.  Teeter  and  Walter 
Pankow,  Estuarine  Processes  Branch,  and  edited  by  Mrs.  M.  C.  Gay,  Information 
Technology  Laboratory,  WES. 

Commander  and  Director  of  WES  during  preparation  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


inches 

knots  (international) 
miles  (US  statute) 
pounds  (force)  per  square  inch 
pounds  (mass)  per  cubic  foot 
tons  (2,000  pounds,  mass) 


By  _ To  Obtain _ 

0.02831685  cubic  metres 

0.3048  metres 

2.540  centimetres 

0.5144444  metres  per  second 
1.609347  kilometres 

703.1  kilograms  per  square  metre 

16.01846  kilograms  per  cubic  metre 

907.1847  kilograms 
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THE  ATCHAFALAYA  RIVER  DELTA 
FIELD  DATA 

Section  2.  Settling  Characteristics  of  Bay  Sediments 

PART  I:  INTRODUCTION 

Background 

1.  The  main  energy  source  to  the  Atchafalaya  Bay  system  is  the 
Atchafalaya  River.  The  river  captures  about  30  percent  of  the  latitude  flow 
(combined  flow  of  the  Mississippi  River  and  Red  River  at  the  latitude  of 

31°  N)  at  the  Old  River  Diversion  Structure  (Figure  1),  and  carries  with  it  an 
average  of  about  100  million  tons*  of  sediment  in  suspension  each  year.  Over 
the  past  several  decades ,  the  suspended  sediment  has  filled  in  the  Atchafalaya 
basin  floodway  between  its  natural  levee  systems  and  is  now  depositing  rapidly 
in  Atchafalaya  Bay  (Figure  1).  As  shown,  two  deltas  are  forming  in 
Atchafalaya  Bay:  at  the  mouth  of  Lower  Atchafalaya  River  and  at  Wax  Lake  Out¬ 
let  (WLO) .  The  evolving  deltas  became  subaerial  and  vegetated  in  1973  and 
have  since  become  one  of  the  most  dynamic  currently  active  delta  systems  in 
the  world.  The  evolving  deltas  have  converted  shallow  bays  into  marshes  and 
continue  to  generate  a  great  deal  of  interest  in  deltaic  processes.  The  pri¬ 
mary  benefit  from  these  two  deltas  has  been  the  addition  of  new  land  to  the 
coast  of  Louisiana  in  areas  otherwise  experiencing  land  loss.  The  primary 
concerns  with  the  evolving  deltas  have  been  siltation  ir.  the  navigation  chan¬ 
nels  and  backwater  flooding  in  the  surrounding  low-lying  coastal  parishes  of 
southern  Louisiana.  Fine-grained  sediments  predominate  in  the  Atchafalaya  Bay 
system. 

2.  Phenomenal  growth  of  the  subaerial  Atchafalaya  River  delta  and  the 
emerging  WLO  delta  led  the  US  Army  Engineer  District,  New  Orleans,  to  request 
that  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  conduct  an  inves¬ 
tigation  to  predict  future  growth  of  the  deltas  and  effects  of  that  growth. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 
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Figure  1.  Vicinity  sketch  showing  the  Atchafalaya  Rive 
and  Wax  Lake  Outlet  deltas 


Objective 


3.  The  overall  objective  of  the  Atchafalaya  Bay  investigation  is  to 
predict  future  evolution  of  the  delta  system.  Specific  answers  are  sought  to 
these  questions: 

a.  For  existing  conditions  and  no  actions  other  than  those  already 
practiced  (i.e.,  maintenance  of  navigation  channels),  how  will 
delta  deposits  evolve  over  the  short-to-medium  term 

(10-15  years)  and  the  long  term  (50  years)? 

b.  How  will  the  delta's  evolution  affect: 

(1)  Flood  stages? 

(2)  Maintenance  dredging  of  the  navigation  channel? 

13)  Salinity,  sedimentation,  and  circulation  in  the 
Atchafalaya  Bay  system? 

c.  What  would  be  the  impact  of  various  navigation/flood-control 
structural  alternatives  on  all  of  the  above? 

4.  This  report  summarizes  studies  conducted  by  the  WES  Hydraulics  Labo¬ 
ratory  to  characterize  sediments  depositing  in  Atchafalaya  Bay,  Louisiana. 
These  studies  were  undertaken  as  part  of  a  comprehensive  investigation  of  the 
deltaic  evolution  occurring  in  Atchafalaya  Bay.  Objectives  were  to  describe 
field  conditions  and  to  generate  information  needed  by  the  WES  numerical  sedi¬ 
ment  transport  modeling  effort  for  the  prediction  of  delta  evolution.  The 
depositional  properties  of  sediments  were  the  primary  subject.  These  proper¬ 
ties  include  the  settling  velocity  distributions  of  sediments  found  in  suspen¬ 
sion  and  in  bed  deposits,  the  concentration  of  newly  deposited  sediments,  and 
the  rates  of  consolidation  of  newly  deposited  bed  sediments. 

Approach 

5.  The  plan  of  the  overall  investigation  includes  multiple  techniques 
to  predict  delta  growth.  Each  method  builds  upon  prior  work  and  employs  pre¬ 
dictive  techniques  using  progressively  greater  degrees  of  sophistication, 
including  the  following: 

a.  Extrapolation  of  observed  bathymetric  changes  into  the  future. 

b.  A  generic  analysis  that  predicts  future  growth  by  constructing 
an  analogy  between  behavior  of  the  Atchafalaya  Delta  and  other 
deltas  in  similar  environments. 
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Sampling  stations  in  the  Atchafalaya  Bay  system,  Louisiana 


properties  of  the  sediments  themselves,  conditions  under  which  deposition  and 
erosion  occur  in  Atchafalaya  Bay  are  of  interest  for  a  general  understanding 
of  sedimentat ion  processes  in  the  system  and  also  for  numerical  sediment 
modeling  efiorts.  Data  on  flows  and  sediment  concentration  distributions  were 
used  to  identify  these  conditions.  Separate  laboratory  studies  were  also  con¬ 
ducted  on  this  subject  (Dixit,  Mehta,  and  Partheniades  1982).  Data  on  the 
grain  size  and  radiograph  analysis  of  sediment  samples  collected  throughout 
the  bay  system  are  presented  in  Report  2,  Section  3,  of  this  report  series  fin 
preparation) . 

9.  Settling  tests  are  covered  in  Part  II  along  with  data  on  the  initial 
bed  densities  of  sediment  settling  from  suspension.  Consolidation  of  bed 
sediments  in  dense  suspensions  by  hindered  settling  consolidation  is  covered 
in  Part  III.  Part  IV  covers  field  data  on  the  distributions  oi  TSX  and  set¬ 
tling  velocity,  vertical  TSM  stratification,  bed  shear  stress,  and  the  bulk 
wet  densities  and  cation  exchange  capacities  of  bed  sediments. 
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PART  II:  SETTLING  TESTS 


10.  Settling  tests  were  performed  in  settling  tubes  and  by  numerical 
methods.  Settling  tube  tests  are  indirect  methods  of  determining  settling 
velocity.  The  settling  characteristics  of  quiescent  suspensions  were  observed 
over  time  in  the  laboratory  and  in  the  field.  Numerical  models  added  turbu¬ 
lence  and  varying  bottom  boundary  conditions  to  the  study  of  the  settling  of 
suspensions.  Because  sediments  were  tested  under  various  concentrations, 
salinities,  heights,  and  turbulence,  the  laboratory  behavior  of  sediments  was 
observed  and  extended  to  the  field  behavior. 

Basic  Concepts  and  Definitions 

11.  Several  interrelated  factors  influence  the  behavior  of  sediment 
particles  suspended  in  a  fluid.  The  factors  can  be  categorized  in  two  groups: 
the  properties  of  the  sediment  and  the  conditions  in  which  the  sediment 
exists.  Both  of  these  groups  are  important  in  the  classification  scheme  for 
suspended  sediment  and  can  be  studied  in  the  laboratory  using  a  quiescent 
settling  tube.  For  some  purposes,  the  sediment  grain  size  in  itself  is  a  use¬ 
ful  classification  scheme;  however,  a  more  useful  classifier  for  suspended 
sediment  is  settling  velocity.  Basically,  settling  velocity  is  caused  by  the 
imbalance  of  gravitational  forces  and  the  viscous  drag  experienced  by  sus¬ 
pended  particles  and  aggregations  of  particles.  Settling  velocity  takes  into 
account  the  effects  of  shape,  density,  and  physical-chemical  characteristics 
in  addition  to  the  particle  size.  To  classify  the  behavior  of  particles  in 
turbulent  suspensions,  certain  properties  of  the  flow  must  also  be  considered, 
namely,  depth  and  turbulence.  The  particle  Peclet  number  combines  the  effects 
of  settling  and  flow  conditions.  This  number  is  the  ratio  of  settling  time  to 
vertical  mixing  time  and  can  be  related  to  suspension  stratification  and  mass 
flux  to  the  bed. 

12.  Particle  grain  size  distribution  and  mineral  composition  are  im¬ 
portant  factors  in  the  settling  velocity  of  fine-grained  sediments.  Other 
conditions  that  affect  the  settling  velocity  of  fine-grained  sediment  include 
concentration,  settling  height,  salinity,  temperature,  and  turbulence.  Most 
of  these  conditions  can  be  reproduced  in  the  laboratory.  After  testing,  set¬ 
tling  velocities  can  be  converted  into  equivalent  quartz  spherical  diameters 
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for  comparison  to  mineral  grain  sizes  by  using  Stokes'  law  of  settling. 


Examples  of  such  equivalent  diameters  for  a  particle  density  of  2.65  g/cu  cm 


are  shown  in  the  following  tabulation: 
Settling 

Velocity,  mm/sec 


Stokes  Equivalent 
Diameter ,  pm 


0.0009  1 
0.0038  2 
0.0160  4 
0.0600  8 
0.2400  16 
0.9800  32 
3.8000  64 


Figure  3  shows  an  example  comparison  of  grain  size  distribution  (determined  by 
hydrometer  test)  plotted  with  a  series  of  settling  tests  made  on  bed  material 
from  station  F  with  varying  concentrations.  The  largest  differences  are  in 
the  slow-settling  or  fine-grained  range.  These  small  particles  have  high 


tzeja/  va/ent  S/oA res'  /Piamc-fcr ,  rn 

/  2  4  3  /<*  32  <24 


Figure  3.  Comparison  between  grain  size  and  settling  velocity 
distributions  for  station  F 


surface-area-to-volurae  ratios  and  surface  electrical  charge  densities.  If 
conditions  bring  them  into  close  contact  or  if  ionic  concentrations  reduce 
repulsive  forces,  they  will  adhere  to  one  another.  The  strength  of  this  ad¬ 
hesion  depends  on  interparticle  distances.  In  the  example  shown  in  Figure  3, 
concentration  has  affected  the  settling  velocity  distribution.  Thus  the  set¬ 
tling  velocity  distribution  cannot  be  deduced  directly  from  dispersed  grain 
size  distribution. 

13.  Dispersed  grain  size  distributions  for  stations  F  and  M  are  shown 
in  Figure  4.  They  are  characteristic  of  the  area.  A  compilation  of  grain 
size  distributions  for  hundreds  of  locations  in  the  bay  system  is  given  in 
Report  2,  Section  3,  of  this  report  series  (Donnell,  in  preparation).  Bed 
sediments  were  found  generally  to  be  silty,  fine-grained  materials  with 
varying  amounts  of  clay  and  sand. 


Figure  4.  Dispersed  grain  size  distribution  for  two  sampling  stations 

14.  To  represent  mathematically  the  vertical  settling  of  fine  sedi¬ 
ments  in  the  presence  of  turbulence,  the  one-dimensional  transport  equation 
with  Fickian  diffusion  and  diffusivity  assumptions  was  used: 
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ac  _  ac  =  a_ 

3T  '  s  3Z  "  ?2 


(1) 


where 

C  =  concentration 
T  =  time 

W  =  terminal  scalar  settling  velocity  resulting  from  the  balance  of 
gravitational  forces  and  viscous  drag  (the  settling  velocity  in 
quiescent  fluid) 

2  =  vertical  dimension 

=  effective  vertical  eddy  diffusivity 

Since  fine  sediments  are  being  considered,  it  is  assumed  that  is  the  same 

as  the  combined  turbulent  and  molecular  dif fusivities  for  momentum.  The  ver¬ 
tical  dimension  (Z)  is  measured  upward  from  the  bed.  The  following  flux 
boundary  conditions  are  imposed: 


and  either 


or 


or 


w  c  + 

s 

ac 

KZ  az  = 

0 

at 
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(water  surface) 

(2) 

ac 
KZ  az 
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z 
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for  free  settling 

(3) 

w  c  + 

s 

k 

z  az 

0 
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Z  =  0 

(bed)  for  no-settling 

(4) 

w  c  + 

s 

K  1C  = 

z  az 

m 

at 

Z  =  0 

(bed)  for  erosion 

(5) 

where  m  is  the  mass  rate  of  erosion  and  H  is  the  total  depth  or  height. 
Equation  1  can  be  nondimensionalized  by  substituting 


z 


Z 

H  ’ 


t 


W  T 
s 

H  ’ 


and  P 

e 


where  Pg  is  the  particle  Peclet  number.  Then, 


ac  ac  i_  a^c 

9t  "  82  =  pe  az2 


(6) 
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is  a  nondimensional  transport  equation  and  the  nondimensional  flux  becomes 


C 


1_  oC 

P  aZ 
e 


This  implies  that  at  steady  state  or  equilibrium 

P  =  I  5C 

e  C  dz 


(7) 


where  8C/3z  is  normally  negative.  The  nondimensional  parameter  governing 

stratification  is,  therefore,  the  particle  Peclet  number.  This  was  confirmed 

by  numerical  simulations  and  the  relationship  of  Equation  7  verified.  For  a 

given  probability  that  a  particle  reaching  the  bottom  stays  there,  P  af- 

e 

fects  the  rate  of  mass  flux  to  the  bed.  Parametric  tests  were  run  on  the 

effects  of  P^  and  bottom  boundary  conditions  on  stratification,  the  effect 

of  P  on  mass  flux  to  the  bed,  and  the  effects  of  P  on  the  interpretation 
e  e 

of  settling  velocity  analyses.  Depth-averaged  and  parabolically  distributed 
values  of  were  tested  and  compared. 

Description  of  Settling  Tests 


15.  Settling  tube  tests  were  performed  on  suspended  samples  taken  in 
the  field.  Parametric  tests  were  also  run  on  suspended  sediments  and  resus¬ 
pended  bottom  sediments  to  establish  the  effects  of  certain  conditions  on  set¬ 
tling  velocity.  Conditions  tested  included  concentration,  settling  height, 
salinity,  and  turbulence. 

Laboratory  methods 

16.  Tests  were  performed  in  two  clear  acrylic  tubes,  each  2  m  high  by 
10  cm  in  diameter,  with  sampling  ports  constructed  in  them.  One  tube  had 
silicon  rubber  ports  or  diaphragms  and  samples  were  drawn  using  an  18-gage 
hypodermic  needle  inserted  through  the  port  into  the  suspension.  The  other 
tube  had  spigots  with  0.32-cm  bore  diameter  extending  1.5  cm  into  the 
suspension. 

17.  Tests  were  carried  out  by  sampling  from  the  settling  tube  with 
time.  The  sediment  to  be  tested  was  shaken  by  hand  in  a  20-1  carboy  for 
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approximately  3  min  before  being  poured  into  the  settling  tube.  A  sample  was 
taken  immediately,  and  further  sampling  was  spaced  approximately  geometrically 
over  a  6-hr  period.  Between  6  and  11  samples  were  taken  during  each  test. 

TSM  was  determined  on  samples  by  filtration.  Nuclepore  polycarbonate  filters 
with  0.40-pm  pore  size  were  used.  They  were  desiccated  and  preweighed.  A 
vacuum  system  (8-psi  vacuum  maximum)  was  used  to  draw  the  sample  through  fil¬ 
ters.  Filters  and  holders  were  then  washed  with  distilled  water  to  remove 
salts.  Filters  were  dried  at  105°  C  for  1  hr  and  reweighed.  TSM  was  then 
calculated  based  on  the  net  weight  retained  on  the  filter  and  the  volume  of 
the  sample  filtered. 

18.  Some  settling  tests  were  made  at  high  concentration .  In  addition 
to  the  samples  normally  taken,  the  height  of  the  sediment  bed  layer  was  also 
measured,  initially  at  intervals  of  about  2  min.  After  the  settling  velocity 
was  computed,  the  density  of  the  bed  deposits  could  be  estimated  as  the  mass 
of  sediment  settled  divided  by  the  volume  it  occupied  at  a  point  in  time  when 
an  arbitrary  68  percent  of  the  sediment  had  settled  to  the  bed  or  when  the 
test  was  terminated. 

Field  methods 

19.  Field  settling  tests  were  performed  much  the  same  as  laboratory 
experiments  except  that  samples  were  drawn  from  a  field  sampler  immediately 
after  the  sample  was  taken.  Field  stations  are  shown  in  Figure  2.  The  field 
sampler  (and  settling  tube)  was  a  Niskin  type  made  of  opaque  polyvinyl  chlo¬ 
ride.  The  sampler  is  an  0.85-m-long  tube  that  opens  at  both  ends  and  is  held 
horizontally,  aligned  with  the  flow,  for  sampling.  The  field  settling  sampler 
is  shown  in  Figure  5.  After  closure,  it  was  brought  to  the  surface  and  placed 
vertically  in  a  rack,  and  sampling  continued.  Normally,  six  samples  were 
taken  over  a  4-hr  period.  The  field  method  minimizes  the  disturbance  that 
suspended  particles  and  aggregates  experience  during  the  settling  test.  The 
estimates  of  settling  velocity  obtained  in  this  way  are  expected  to  be  closer 
to  in  situ  values  than  those  made  using  laboratory  tests. 

Settling  data  reduction 

20.  Several  methods  were  used  to  analyze  the  settling  tube  data,  de¬ 
pending  on  the  purpose  of  the  tests.  The  accumulation  method  was  applied  to 
multiple-depth  tests.  The  accumulation  method  was  developed  by  Oden  (1925)  in 
1925  and  recommended  for  settling  velocity  determinations  at  low  concentra¬ 
tions  (Inter-Agency  River  Basin  Committee  1953).  This  method  relates  settling 
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Figure  5.  Settling  velocity  field  sampling  tube 

velocity  to  the  differentiation  of  the  sedimentation  curve  with  respect  to 
time.  Usually  a  bottom  withdrawal  tube  is  used  in  tests.  However,  for  this 
study,  the  sedimentation  curves  were  constructed  by  integrating  the  concentra 
tions  from  several  heights  in  a  settling  tube.  Results  from  the  multiple- 
depth  tests  were  also  used  to  construct  concentration  profile  diagrams  using 
method  proposed  by  McLaughlin  (1959).  These  diagrams  can  be  used  to  detect 
the  effects  of  flocculation,  aggregation,  or  hindered  settling  by  examining 
the  slope  of  the  isopleths  of  settling  height  H  divided  by  time  T  .  Most 
routine  settling  velocity  determinations  used  the  pipette  analyses,  assuming 
that  the  change  of  concentration  at  H  distance  below  the  free  surface  repre 
sents  that  fraction  of  the  suspension  with  settling  velocity  greater  than 
H/T  . 

21.  A  computer  program  named  SETL  was  developed  to  analyze  settling 
test  results.  Data  from  the  settling  tube  tests  were  least-squares-f it  by  a 
second-order  polynomial  of  the  form: 

%  settled  =  a  +  a, (In  T)  +  a0  (In  T)^  (8) 

o  i  l 

where  the  a  terms  are  regression  constants.  For  the  accumulation  method, 
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settling  velocity  distributions  were  calculated  by 


%  greater  than  H/T  =  (aQ  -  a^)  +  (a^  +  2&q)  In  T  +  a^Cln  T)2  (9) 

This  technique  is  equivalent  to  older  graphical  methods. 

22.  A  comparison  between  the  regression-method  computer  program  SETT 
and  graphical  techniques  is  shown  in  Figure  6.  The  data  and  graphical  results 
are  from  Christodoulou,  Leimleuhier,  and  Ippen  (1974).  For  the  pipette  anal¬ 
yses,  the  percent  greater  than  H/T  was  computed  direct1 v  from  Equation  8, 
while  increments  of  the  distribution  were  found  by  quadratic  solution  of  that 
equation. 

23.  In  addition  to  the  cumulative  distribution  of  settling  velocity, 
some  statistical  parameters  were  calculated.  They  include  geometric  mean, 
which  is  more  representative  of  the  distribution  than  the  median,  and  the 
dimensionless  geometric  standard  deviation  and  skewness.  The  equations  used 
are  those  suggested  by  Inman  (1963)  for  the  generally  log-normal  particle  size 
distributions : 


Geometric  mean 


W  \ 
16  S84/ 


1/2 


Geometric  standard  deviation 


(10) 


(ID 


Skewness 


log  a 

g 


(12) 


where  W  is  the  settling  velocity  at  the  x  percentile.  Differential 
settling  classes  were  defined  by  a  geometric  series  constructed  by  raising  an 
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a.  Computed  kaolinite  and  field  data 


b.  Boston  Harbor  mud  and  field  data 

Figure  6.  Comparison  between  regression  program 
SETL  and  graphical  method  (MIT  data  obtained  from 
Christodoulou ,  Leimleuhler,  and  Ippen  1974) 
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arbitrary  number  to  progressively  higher  powers: 


(3  *  10  ^  m/sec)^  ,  k  =  1  to  8 

and  the  concentrations  of  suspended  material  falling  within  these  classes  were 
calculated  by  the  computer  program  SETL.  Settling  velocities  were  also  calcu¬ 
lated  at  10  percent  intervals  over  the  settling  velocity  distribution. 

24.  Several  comparisons  among  the  various  settling  analyses  were  made 
in  addition  to  the  comparison  between  graphical  and  regression  model  accumu¬ 
lation  methods  already  mentioned.  Figure  7  shows  an  example  plot  of  a 
multiple-depth  sedimentation  test.  The  data  were  integrated  and  analyzed  by 
the  accumulation  method  as  well  as  individually  analyzed  by  the  pipette 
method.  Results  are  compared  in  the  following  tabulation: 

Fiercer?/  removed  (by  weiyhf') 


_ Settling  Velocity,  mm/sec _ 

Pipette _  _ Accumulation 


Depth,  m 

Median 

Mean 

Standard 

Deviation 

Median 

Mean 

Standard 

Deviation 

0.20 

0.011 

0.015 

5.56 

0.95 

0.010 

0.016 

13.04 

1.65 

0.012 

0.018 

9.47 

0.012 

0.018 

7.66 
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Results  showed  that  the  accumulation  and  pipette  methods  of  settling  velocity 
analysis  were  equivalent,  and  that  settling  velocity  results  from  this  study 
should  be  comparable  to  other  reliable  results. 

25.  According  to  McLaughlin's  interpretation  of  multiple-depth  tests 
using  settling  depth  versus  time  (Z/T)  (McLaughlin  1959),  the  isopleths 
sloping  toward  the  right  with  increased  depth  indicate  an  increase  in  settling 
velocity  with  depth  due  to  aggregation  or  flocculation.  Figure  7  therefore 
indicates  a  slight  increase  in  Wg  with  depth.  The  same  indication  is  given 
in  the  preceding  tabulation  by  the  pipette  analyses. 

26.  The  pipette  analysis  and  computer  program  SETL  were  also  checked  by 

analyzing  numerical  settling  experiments  at  high  numbers.  For  a 

specified  settling  velocity  of  0.200  mm/sec,  analysis  of  numerical  results  by 
the  pipette  method  and  computer  program  SETL  recovered  a  mean  settling  veloc¬ 
ity  of  0.205  mm/sec  and  geometric  standard  deviation  of  1.2  at  =  100  . 

This  test  was  a  further  indication  that  the  pipette  metnod  and  computer  pro¬ 
gram  SETL  accurately  determine  settling  velocity. 

Field/laboratory  test  difference 

27.  Figure  8  shows  a  comparison  between  field  and  laboratory  settling 
tests  run  on  samples  taken  at  the  same  location  with  the  same  sampler  and  as 
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Figure  8.  Comparison  between  field  and  laboratory 
settling  tests 
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close  in  time  as  possible.  That  comparison  shows  that  the  field  results  were 
higher  than  those  obtained  with  the  laboratory  settling  tube.  The  laboratory 
tests  included  two  different  settling  heights.  Similar  results  have  been  re¬ 
ported  by  other  investigators  (Owen  1971).  The  explanation  given  was  that 
natural  turbulence  brings  aggregates  to  a  state  at  which  they  will  settle  more 
quickly.  As  will  be  discussed  later,  settling  height  also  increases  settling 
velocities.  When  a  small  portion  of  a  relatively  deep  suspension  is  col¬ 
lected,  the  field  settling  tests  are  affected  by  the  in  situ  aggregation 
caused  by  both  turbulence  and  settling  height. 

Numerical  methods 

28.  The  governing  equation.  Equation  1,  along  with  appropriate  boundary 
and  initial  conditions,  was  solved  using  a  method  of  weighted  residuals. 
Specifically,  the  method  of  orthogonal  collocation  was  used  where  the  weight¬ 
ing  function  is  a  Dirac  delta  function  and  the  trial  functions  are  a  series  of 
Legendre  polynomials.  The  collocation  points  or  nodes  are  located  at  Gauss 
points.  Instead  of  breaking  the  domain  into  subelements,  a  single  high-order 
element  with  seven  nodes  was  used.  This  was  easy  to  apply  and  maintained  con¬ 
tinuous  first  and  second  derivatives  throughout  the  domain.  A  fourth-order 
Runge-FCutta  method  was  used  to  integrate  in  time.  Gaussian  quadrature  was 
used  to  integrate  over  the  spatial  domain.  The  interpolated  values  used  to 
construct  and  plot  concentration  profiles  were  generated  from  nodal  values 
using  the  trial  functions  for  interpolation. 

Test  Results 


29.  Table  2  is  a  tabulation  of  the  field  and  laboratory  settling  tests 
performed  by  sampling  station.  The  table  also  shows  other  analyses  performed, 
which  will  be  more  thoroughly  discussed  in  Part  IV.  Settling  velocity  distri¬ 
butions  for  field  settling  tests  are  given  in  Table  3.  Table  4  gives  a  sum¬ 
mary  of  the  turbulence  settling  test  conditions  and  results.  Table  5  summa¬ 
rizes  the  in^ial  bed  densities  observed  in  high-concentration  settling 
experiments. 

Discussion  of  laboratory  Results 


30.  The  field  settling  tests  will  be  discussed  in  Part  IV  in  the 
context  of  other  field  measurements. 
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Effect  of  concentration 

31.  Increased  concentration  decreases  interparticle  distances  and  gen¬ 
erally  enhances  aggregation.  There  is  an  upper  limit  to  this  effect  above 
which  settling  is  hindered.  This  will  be  further  discussed  in  Part  III.  Fig¬ 
ure  9  shows  the  effect  on  bed  samples  from  stations  H  and  F.  Settling  veloci¬ 
ties  are  plotted  for  25,  50,  and  75  percentiles.  Figure  10  shows  the  effect 
of  concentration  on  the  settling  velocity  of  bed  sediments  from  station  M. 

The  greatest  effect  of  concentration  on  Wg  was  found  in  the  75-percentile 
fraction  of  station  H.  Station  F  had  somewhat  less,  and  station  M  showed  lit¬ 
tle  dependence  of  settling  velocity  at  low  concentrations  (below  1  g/£) ,  and  a 
moderate  dependence  between  1  and  10  g /£.  Cation  exchange  capacities  (CEC)  of 
samples  from  stations  H,  F,  and  M  were  32,  27,  11  meq/100  g,  respectively. 
Therefore,  CEC  correlated  with  the  concentration  effects  found  in  the  settling 
experiments.  For  these  samples,  the  higher  the  CEC,  the  greater  the  effect  of 
concentration  on  settling  velocity. 

Effects  of  salinity 

32.  Suspended  samples  were  tested  with  increasing  salinity.  Increased 
ionic  concentration  reduces  electrorepulsive  charges  around  the  particles, 
allowing  aggregation  to  proceed.  Figure  11  shows  the  results  of  a  test  that 
consisted  of  a  mixture  of  suspended  samples  from  stations  C  and  M.  Native 
ocean  water  was  used  to  increase  the  salinity.  Suspended  samples  from  sta¬ 
tions  H  and  F  were  tested  and  results  shown  in  Figure  12.  Their  salinity  was 
increased  using  a  saturated  brine  solution  made  from  Louisiana  salt  dome  rock 
salt . 

33.  Composite  sediments  from  stations  C  and  M  showed  an  increase  in  the 
median  settling  velocity  by  a  factor  of  about  2.5  when  the  salinity  was  in¬ 
creased  from  about  0.25  ppt  (river  concentration)  to  between  2.3  and  2.2  ppt. 
Station  H  showed  little  change  while  station  F  actually  showed  a  decrease  in 
settling  velocity  between  0.25  and  2  ppt  followed  by  a  large  increase  in 
settling  velocity  at  over  30  ppt.  Experimental  error  is  suspected  in  the 
2-ppt  test. 

Effects  of  settling  height 

34.  Some  tests  on  the  effects  of  settling  height  on  settling  velocity 
were  made.  Differential  settling  can  lead  to  further  aggregation  and  in¬ 
creased  settling  velocity.  Figure  13,  which  shows  test  results  for  stations  H 
and  M  at  different  settling  heights,  indicates  increased  settling  velocity 
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m/  sec 


mm  / 5  <?c. 


Not  a  :  River  sa/j'ni+y  assumed  a  /  0. 2£ pp  f 

Figure  11.  Effects  of  salinity  on  the  settling  of 
stations  C  and  M  sediments 


Figure  12.  Effects  of  salinity  on  the  settling  of 
stations  F  and  H  sediments 


Figure  13.  Effects  of  settling  height  on 
settling  velocity 

with  settling  height.  It  is  generally  reported  that  this  effect  ceases  over 
about  a  2.0-m  height.  This  was  not  tested.  Laboratory  vertical  clearance 
limited  settling  tests  to  1.8-m  settling  heights.  As  mentioned,  field  tests 
used  0.85-m  tubes. 

Effects  of  turbulence 

35.  Turbulence  can  have  different  effects  on  sediment  aggregates  and 
their  settling  velocities.  Turbulent  shears  can  drive  particles  closer  to¬ 
gether,  overcoming  repulsive  forces,  and  thus  form  relatively  dense,  strong, 
faf-jr  settling  aggregates.  Turbulence  can  also  destroy  aggregates  that  are 
too  weak  to  withstand  the  imposed  shear.  Tests  were  performed  to  determine 
which  of  these  effects  dominate  under  natural  conditions. 

36.  Laboratory  settling  tube  tests  were  performed  with  an  oscillating 
screen  used  to  generate  turbulence.  A  screen  of  2-mm  plastic  mesh  was  at¬ 
tached  by  a  long  rod  in  the  100-mm-diameter  tube  and  attached  to  a  variable- 
speed  oscillator  with  a  19-mm  stroke  mounted  above  the  settling  tube.  The 
screen  was  centered  between  the  free  surface  and  the  sampling  point.  A  1-m 
settling  height  was  used  for  these  tests.  The  sampling  port  was  located  25  cm 
above  the  physical  bottom  boundary  of  the  tube.  Three  oscillation  speeds  were 
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used  and  measured  using  a  touchless  tachometer. 

37.  Flow  visualization  tests  made  with  the  screen  oscillating  indicated 
that  turbulence  was  most  intense  at  the  screen  and  decayed  away  from  it.  At 
the  highest  oscillating  speed,  turbulence  was  noticeable  throughout  the  set¬ 
tling  depth  but  decayed  to  very  low  levels  near  the  physical  bed  of  the  tube. 
Resuspension  of  material  that  settled  to  the  bottom  was,  therefore,  unlikely. 

38.  Dye  experiments  were  performed  in  the  settling  tube  and  the 

2 

change  of  dye  variance  (a  )  with  time  used  to  estimate  dif fusivities 

2 

=  l/2(3a  /9T) .  For  the  sizes  of  material  tested,  the  analogy  between  dif- 
fusivities  of  conservative  dye  and  sediment  was  valid.  Depth-averaged  dif- 
fusivities  estimated  are  shown  in  the  following  tabulation: 


Screen 
Oscillation 
Rate,  rpm 

90 

250 

440 


Dif fusivity 
2  . 

mm  /sec 
45 
160 
285 


Except  for  the  oscillating  screen,  settling  tests  were  otherwise  run  normally. 

39.  As  mentioned  earlier,  settling  tube  tests  imply  settling  velocities 

from  the  sedimentation  characteristics,  specifically  the  time-history  of  con¬ 
centration.  Results  of  the  turbulence  tests,  as  shown  in  Table  4,  indicated 
that  sedimentation  characteristics  of  the  suspensions  were  generally  changed. 
Effective  values  of  decreased  with  increased  oscillation  speeds.  This 

does  not  mean  that  the  settling  velocity  of  the  particles  themselves  changed, 
however . 

40.  To  clarify  the  effect  of  turbulence  on  the  particle  aggregates, 
numerical  settling  experiments  were  performed.  Using  this  method,  the  effect 
of  turbulence  on  sedimentation  was  separated  from  effects  on  individual  aggre¬ 
gates.  The  numerical  methods  were  described  earlier.  Numerical  settling  ex¬ 
periments  provided  time-histories  of  suspension  concentration  for  a  fixed 

W  .  The  same  analyses  were  performed  on  the  numerical  results  as  were  done 
on  the  experimental  values  to  estimate  mean  effective  settling  velocities. 
Results  were  used  to  plot  mean  effective  settling  velocities  normalized  to  the 
90-rpm  value  versus  a  particle  Pg  calculated  using  the  90-rpm  Wg  and 
for  the  test.  If  the  actual  Wg  were  unchanged  with  increased  turbulence, 
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then  the  effective  Wg  measured  in  the  laboratory  should  be  reduced  in  the 

same  manner  as  the  numerical  results,  where  W  was  held  constant. 

s 

41.  The  numerical  results  are  plotted  in  Figure  14  by  normalizing  the 
effective  W  to  P  =  1  .  Since  the  values  of  P  for  Test  OE  were  higher 

S  6  6 

than  for  other  tests,  numerical  results  are  also  plotted  unnormalized  for  com¬ 
parison  to  this  test.  Although  there  is  scatter,  the  slope  of  the  numerical 
results  is  close  to  Tests  OB,  OA,  and  OE.  Tests  OF  and  OD  had  smaller  slopes, 
and  this  may  indicate  that  mild  turbulence  had  a  positive  effect  on  . 
Overall,  however,  the  range  of  turbulence  used  in  these  tests  had  no  dramatic 
effect  on  the  particle  or  aggregate  settling  velocities.  The  experimental  P^ 
range  overlaps  the  field  measured  P^  ,  but  lower  values  in  the  field-measured 
P^  indicate  more  turbulence  might  be  applied  to  future  laboratory  studies. 
Initial  bed  density 

42.  Concentration  and  densities  of  initial  bed  deposits  from  resus¬ 
pended  station  M  sediments  were  similar  for  the  range  covered  by  the  high- 

concentration  settling  tests.  As  noted  previously,  at  concentrations  above 

3 

the  settling  range  (>.130  g/2.  or  0.13  g/cm  dry  weight  for  station  M)  ,  hindered 
settling  consolidation  occurs  and  the  sediment  settles  or  consolidates  in 
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Figure  14. 


Effects  of  turbulent  P^  on  effective  settling 
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mass.  The  experiments  described  here  are  stiil  in  the  settling  range  of  con¬ 
centration  so  that  as  settling  proceeds,  a  distinct  dense  layer  forms  below 

the  suspension  at  the  bottom  of  the  settling  tube.  Composite  sediments  from 

3 

station  M  were  formed  with  initial  bed  densities  of  about  1.45  g/cm  bulk  wet 

3  3 

density  (BWD)  or  0.72  g/cm  dry  weight  compared  with  1.8  g/cm  BWD  or 

3 

1.27  g/cm  dry  weight  measured  in  field  samples.  Field  measurements  are 

presented  in  Part  IV.  The  fine  fraction  of  station  M  sediments  had  much  lower 

initial  bed  densities  than  would  be  expected.  This  material  is  more  cohesive 

and  characteristically  would  be  expected  to  have  low  initial  densities.  Still 

3  3 

they  averaged  about  1.24  g/cm  BWD  or  0.39  g/cm  dry  weight  including  one  low 
value . 

Numerical  stratification  tests 

43.  Stratification  was  evaluated  as  6C/C  where  6C  is  the  surface- 

to-bottom  difference  in  concentration  and  C  is  the  depth-averaged  concentra¬ 
tion.  Model  results  with  free-settling  bottom  boundary  condition  and  depth- 
averaged  compared  favorably  to  published  results  that  used  other 

numerical  methods.  Model  results  showed  that,  as  indicated  by  Equation  7, 
stratification  is  related  to  Pg  for  a  given  bottom  boundary  condition.  For 
a  no-settling  bottom  boundary  condition  and  depth-averaged  ,  6C/C  is 

equal  to  Pg  .  For  the  same  no-settling  boundary  condition  and  parabolically 
distributed  K z  ,  6C/C  =  1.5Pe  over  the  range  0.01  <  <  1  (where  the  par¬ 

ticle  Pfi  is  calculated  based  on  a  depth-averaged  K^).  Since  the  parabolic 
distribution  is  more  realistic  than  a  uniform  distribution,  these  results  were 
used.  For  free-settling  bottom  boundary  condition,  6 C/C  =  0. 75Pe  for 

0.01  <  Pg  <  1  .  Above  P  *  1  ,  these  relationships  diverge  slightly.  Fig¬ 
ure  15  shows  the  numerically  generated  suspension  classification  diagram. 
Figure  16  shows  example  equilibrium  profiles  with  no-settling  bottom  boundary 
conditions  and  Figure  17  shows  free-settling  results  at  Pg  =  0.5  . 

44.  These  results  suggest  that  the  stratification  of  a  suspension 
depends  on  P^  and  bottom  boundary  condition.  Therefore,  if  stratification 
and  particle  P^  can  be  accurately  determined,  erosion  or  deposition  that 
occurs  over  spatial  scales  of  greater  than  perhaps  500  ra  could  be  estimated. 
Comparison  to  field  data  will  be  made  later  in  Part  IV. 

Mass  flux  to  the  bed 

45.  Mass  flux  to  the  bed,  or  deposition  D  ,  can  be  described  as 
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P€  =  Pec/e  t  number 


Figure  15.  Suspension  classif ication  diagram 
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D  =  FW  C 
s 

where  P  is  the  probability  of  a  particle  reaching  the  bed  remaining  there, 

and  C  refers  to  the  concentration  just  above  the  bed.  P  has  been  found  to 

vary  with  bed  shear  stress  t,  ,  and  be  0  at  some  critical  value  t  ,  and 

b  cd 

l  at  x,  =  0  . 
b 

46.  Analytic  solutions  to  sediment  transport  equations,  such  as 
Equation  1,  for  deposition  and  under  certain  simplifying  assumptions  can  be 
made.  Two  such  solutions  were  tested  using  numerical  model  results,  and  the 
effect  of  Pg  on  deposition  evaluated.  The  two  analytic  solutions,  assuming 
Lhac  P  =  1  ,  are 


C 

—  =  eXp 
o 


(14) 


where  C  is  the  initial  concentration  at  T  =  0  . 
o 

47.  Plots  of  these  functions  and  of  numerical  results  at  P  100  and 

e 

P^  =  1  are  shown  in  Figure  18.  Tne  numerical  results  for  P^  =  1  were  ob¬ 
tained  using  a  uniform  initial  concentration.  For  =  100  ,  the  model  rur, 

was  made  with  an  initial  concentration  varying  linearly  from  0  at  the  surface 
to  200  percent  of  the  average  at  the  bottom,  a  more  realistic  condition  for 
this  very  high  P^  . 

48.  Results  showed  that  at  values  of  P^  of  about  1  or  less,  the  expo¬ 

nential  function  described  deposition  quite  well.  As  will  be  shown  later, 
field  were  generally  in  this  range. 
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PART  III:  CONSOLIDATION  TESTS 


49.  After  sediments  first  settle  to  the  bed,  they  undergo  further  set¬ 
tling.  This  is  termed  hindered  settling  consolidation  because  concentrations 
are  in  the  hindered  settling  range  and  thus  settling  rates  decrease  with  in¬ 
creased  concentration.  This  is  significantly  different  from  the  settling 
processes  described  in  Part  II. 

50.  Self-weight  consolidation  and  gravity  compaction  were  not  addressed 
in  this  report.  They  are  important  to  consolidation  after  hindered  settling 
consolidation  has  taken  place,  especially  in  thick  sediment  layers.  Important 
changes  in  density  and  void  ratios  usually  occur  with  the  first  few  hundred 
metres  of  burial.  However,  the  deposits  that  have  accumulated  in  Atchafalaya 
Bay  over  recent  decades  are  on  the  order  of  2-3  m  thick  or  less,  and  in  situ 
densities  have  been  measured. 

51.  In  estuaries,  and  Atchafalaya  Bay  in  particular,  hydrodynamic  con¬ 
ditions  are  quasi-periodic  cr  transient.  Slack-water  periods  with  low  shear 
stresses  on  the  bed  can  allow  sediments  to  deposit.  This  is  followed  '-y 
periods  of  higher  shear  stress  and  possibly  erosion.  The  critical  shear 
stress  for  erosion  of  fine-grained  sediments  has  been  found  to  be  strongly 
dependent  on  density.  To  predict  this  behavior,  it  therefore  becomes  impor¬ 
tant  to  estimate  initial  densities  of  the  settled  material  and  the  increase  in 
density  of  this  material  over  periods  of  hours  to  days.  This  was  the  purpose 
of  the  consolidation  tests. 

Basic  Concepts  and  Definitions 

52.  Hindered  settling  occurs  when  the  concentration  of  a  suspension 
becomes  so  large  that  interparticle  distances  impede  settling.  Hindered  set¬ 
tling  consolidation  occurs  at  still  higher  concentrations  when  particles  be¬ 
come  immobilized  in  the  suspension  and  the  suspension  settles  as  a  mass. 
Hindered  settling  has  been  empirically  described  by 

W  =  W  (1  -  KC)5  (lb) 

h  s 

where 

=  hindered  settling  velocity 
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K  =  an  empirical  coefficient,  g/cm  /dry  weight 
An  exponent  of  4.67  is  sometimes  used. 

53.  To  verify  this  relationship  and  to  obtain  further  insights  on  hin¬ 
dered  settling  consolidation  under  field  conditions,  a  numerical  model  was 
constructed.  The  transport  equation  without  diffusion  was  used. 

H  -  h  <whc)  ■  0  <>7> 

Using  a  nondimensional  depth  term  z  =  Z/H  ,  the  transport  equation  becomes 

i  -  s  h  (whc)  ■ 0  n8) 


where  H  is  the  thickness  or  depth  of  the  consolidating  mud.  The  stationary 
lower  boundary  is  z  =  0.0  ,  and  z  =  1.0  is  the  moving  upper  boundary.  At 
the  lower  boundary  (z  =  0.0)  the  no-flux  condition  W^C  =  0  was  imposed. 

At  the  upper  boundary,  no  boundary  condition  was  imposed  and  the  upper  bound¬ 
ary  was  allowed  to  settle.  The  thickness  of  the  mud  was  calculated  by 


1  .  U 

I  C  dz 

z=0.0 


(19) 


where  m  is  the  mass  per  unit  surface  area  of  the  consolidating  layer,  a 
constant  for  a  given  problem. 

54.  A  solution  method  similar  to  that  presented  in  Part  II  was  used  for 

the  numerical  experiments.  At  each  time-step,  the  variable  coefficients  W, 

h 

and  H  were  calculated  by  Equations  18  and  19.  Normally  an  arbitrary  301 
time-steps  were  made.  Every  tenth  was  printed  and  a  data  file  created.  Plots 
were  made  of  vertical  profiles  of  concentration  with  nondimensional  depth. 
Plots  of  thickness  H  were  made  versus  time. 

Description  of  Tests 

55.  Hindered  settling  consolidation  was  observed  in  the  laboratory  by 
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mixing  bottom  sediments  with  native  water  to  make  a  danse  suspension  or  mud, 
then  observing  its  consolidation  in  a  clear  tube.  Tubes  200  cm  high  by  10  cm 
in  diameter  of  clear  acrylic  and  45  cm  high  by  22  cm  in  diameter  of  clear 
plastic,  both  equipped  with  0.32-cm-bore  Teflon  spigots,  were  used.  At  high 
concentrations,  suspensions  settle  as  a  mass  and  a  very  distinct  clear  layer 
forms  above  the  consolidating  sediment.  Observations  were  made  of  the  loca¬ 
tion  of  this  interface  with  time.  The  rate  of  change  in  the  interface  height 
at  the  beginning  of  the  tests  was  equal  to  .  Samples  were  drawn  with 

depth  through  spigots  and  analyzed  for  BWD  by  pycnometer  method  during  some 
tests. 

56.  BWD  is  the  density  of  sediment  and  water  mixture.  BWD  analyses  are 

3 

performed  using  wide-mouth,  25-cm  ,  constant-volume  pycnometers  calibrated  for 
tare  weight  and  volume.  To  make  a  measurement,  a  pycnometer  was  partially 
filled  with  sediment  and  weighed.  It  was  then  topped  off  with  distilled 
water.  Care  was  taken  to  remove  any  bubbles  from  the  pycnometer  before  it  was 
reweighed.  The  BWD  of  the  sediment  was  then  calculated: 


BWD  = 


(P^Msed  wt  -  tare  wt) 


(p  ) (vol  pyc)  +  (sed  wt)  -  (sed  +  water  wt) 
w 


(20) 


where 


P  = 


w 

sed  wt 
tare  wt 
vol  pyc 
sed  +  water  wt 


density  of  water  at  temperature  of  analysis 
weight  of  pycnometer  with  sediment 
tare  weight  of  pycnometer 
volume  of  pycnometer 

weight  of  pycnometer  st-.liu....,  o.  d  water 


Test  Results  and  Discussion 


57.  Table  6  lists  the  test  material  sample  stations,  initial  condi¬ 
tions,  duration  of  consolidation  experiments,  and  final  average  concentration 
for  tests  in  which  samples  were  drawn.  Plots  of  these  consolidation  tests  are 
given  in  Plates  1-6.  There  was  no  indication  of  a  delay  period  before  hin¬ 
dered  settling  consolidation  began,  as  has  been  reported  for  some  cohesive 
sediments . 

58.  Other  tests  were  performed  in  which  only  the  initial  concentration 
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and  interface  heights  were  measured.  Station  M  sediments  were  used  for  these 
tests,  and  results  are  shown  in  Figure  10.  By  associating  with  C  ,  the 

effect  of  Cq  on  settling  was  determined  for  this  sediment. 

59.  Numerical  model  results  indicate  that  the  ranges  of  the  constants 
and  K  of  0.00005  to  0.00020  m/sec  and  2.4  to  2.9  cm’Vg,  respectively, 

can  be  used  to  represent  the  hindered  settling  consolidation  observed  in  the 
laboratory.  Values  of  W  are  larger  and  K  lower  near  the  Atchafalava  Bay 
river  mouth  where  material  is  coarser  and  less  cohesive. 

60.  Vertical  concentration  profiles  and  interface  locations  for  a  sta¬ 
tion  F  sediment  test  are  plotted  in  Figure  19,  along  with  numerical  model  re¬ 
sults.  The  agreement  in  the  interface  location  is  excellent.  The  near-bottom 
BWD's  were  somewhat  lower  in  the  numerical  model,  but  this  was  probably  due  to 
a  coarse  fraction  in  the  prototype  sediment  that  may  have  been  capable  of 
settling  through  finer  sediments.  The  range  of  concentrations  tested  was  not 
very  far  into  the  hindered  settling  consolidation  range  so  that  it  is  possible 
that  some  particles  were  not  immobilized  in  the  suspension  at  the  onset  ol  the 
test . 

61.  The  inflection  (curvature)  in  the  plot  of  interface  ievel  was  re¬ 
produced  in  the  numerical  model.  Tt  was  thought  that  this  feature  signalled 
the  onset  of  self-weight  consolidation  or  compaction.  The  numerical  model 
showed  this  to  be  the  point  when  the  concentration  or  density  inflection 
reached  the  sediment  surface.  The  consolidation  process  can  be  seen  in  Fig¬ 
ure  19  as  a  propagation  of  densif ication  that  began  at  the  bottom  and  propa¬ 
gated  upward.  Local  gradients  in  W,  and  C  caused  concentrations  to  in- 

h 

crease.  At  the  beginning  of  the  test®  rho  suspension  was  uniform  and  the  only 
gradients  existed  at  the  fixed-bottom  boundary.  Laboratory  tests  also  indi¬ 
cate  that  densif icat ion  began  at  the  bottom  and  propagated  upward.  The  thin¬ 
ner  the  layer,  the  more  quickly  the  surface  density  increased.  Near-surface 
density  increase  corresponded  to  the  inflection  in  the  interface  descent.  For 
a  given  K  and  ,  then,  the  increase  in  concentration  at  a  point  depends 

on  its  distance  above  the  bottom  stationary  level,  not  on  depth  from  the 
interface . 

62.  Interface  heights  observed  in  the  laboratory  and  numerically  cal¬ 
culated  for  a  test  using  sediment  from  station  M  are  shown  in  Figure  20.  The 
agreement  is  good. 

63.  These  results  apply  only  to  laboratory  experiments.  Field 
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Figure  20.  Experimental  and  numerical  consolidation  results 
for  station  M  (Atchafalaya  mud) 


measurements  are  reported  in  Part  IV.  Experimental  results  presented  in 
Part  II  indicate  higher  initial  and  final  bed  densities  than  those  found  in 
the  laboratory  tests.  This  is  due  to  the  difference  in  characteristics  be¬ 
tween  deposits  formed  as  a  result  of  settling  and  those  found  by  the  hindered 
settling  consolidation  of  dense  suspensions.  However,  the  hindered  settling 
equation  applies  to  both  cases,  and  numerical  simulations  were  made  to 
represent  field  conditions. 


PART  IV:  ANALYSIS  OF  FIELD  SEDIMENTATION  PROCESSES 


64-  The  overall  Atchafalaya  Bay  project  conducted  by  WES,  of  which  this 
work  is  a  part,  involved  a  large  data  gathering  effort.  Physical  conditions 
including  tide,  current  speed  and  direction,  wind,  and  wave  height  were  mea¬ 
sured  during  a  2-year  period.  Complete  results  are  given  in  Coleman  et  al. 
(1988),  but  some  are  summarized  here  in  relation  to  suspended  sediment  char¬ 
acteristics.  Field  measurements  of  the  settling  velocity  of  suspended  mate¬ 
rial  were  made  using  methods  already  described  in  Part  II.  Samples  were  also 
collected  for  TSM  analysis.  Often  they  were  collected  along  with  salinity  and 
current  measurements  in  vertical  profiles  or  at  select  locations  in  the  verti¬ 
cal.  Field  equipment  and  procedures  are  given  in  Coleman  et  al.  (1988).  CEC 
and  BWD  were  measured  on  subsamples  of  sediment  cores.  All  these  data  were 
analyzed  to  define  and  describe  field  processes  affecting  erosion,  transport, 
and  deposition  of  sediments. 


Basic  Concepts  and  Definitions 


65.  Vertical  profiles  of  TSM  concentration  and  velocity  can  be  used  to 
calculate  unit  sediment  discharge,  sediment  velocity,  shear  velocity,  sediment 
suspension  stratification,  depth-averaged  sediment  concentration,  and  depth- 
averaged  velocity.  Measurements  of  velocity  taken  in  the  vertical  were  fit  to 
a  regression  equation.  The  following  universal  logarithmic  profile  was  used: 


(21) 


where 

If  =  current  velocity 

1/2 

U*  =  bed  shear  velocity,  (t/p)  where  p  is  the  mass  density  of 

the  fluid 

k  =  assumed  to  be  equal  to  0.4 
Zq  =  roughness  height 

The  regression  was  used  to  estimate  the  shear  velocity,  designated  U* ( f  1 )  , 
and  to  estimate  current  velocity  at  a  number  of  points  (Gauss  points)  in  the 


vertical.  From  the  calculated  depth-averaged  velocity.  Manning's  equation  was 
also  used  to  estimate  the  shear  velocity, 


U*(fm)  =  g1/2nU  H  1/6 


(22) 


where 

U*(fm)  =  shear  velocity  estimated  using  Manning's  equation 
g  =  acceleration  due  to  gravity 
n  =  Manning's  coefficient  taken  as  0.020 
U  =  depth-averaged  current  velocity 

Both  methods  of  estimating  shear  velocity  are  subject  to  uncertainty.  The 
logarithmic  profile  fit  is  sensitive  to  uncertainties  in  the  near-bottom 
velocity  values  and  value  of  .  More  near-bottom  measurement  points  would 

have  reduced  uncertainties,  but  could  not  be  obtained.  Manning's  equation  is 
sensitive  to  its  assumed  coefficient.  For  many  profiles,  however,  reasonable 
agreement  was  found  between  the  two  methods.  Where  the  logarithmic  fit  failed 
completely,  giving  a  widely  different  or  sometimes  even  negative  shear  veloc¬ 
ity,  a  synthetic  data  point  was  added  to  the  profile  at  the  top  of  the  laminar 
sublayer  based  on  the  Manning's  estimated  shear  velocity.  The  depth-averaged 
values  of  velocity  and  concentration  were  calculated  by 

H  H 

U  =  4  /  U  dZ  and  C  =  i  /  C  dZ  (23) 

H  Z=0  H  Z=0 

using  Gaussian  quadrature.  The  discharge  of  sediment  per  unit  width  of  flow 

is 

Q(sed)  =  HC  U(sed)  =  f  CU  dZ  (24) 

Z=0 


where  U(sed)  is  the  concentration-weighted  sediment  velocity 

H 

J  CU  dZ 

U(sed)  =  -  (23) 

f  C  dZ 

z=o 
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The  nondimensional  stratification  (<$C/C)  was  defined  in  Part  II  as  the 
surface-to-bottom  difference  normalized  by  the  depth-averaged  concentration. 

To  estimate  this  parameter,  an  exponential  curve  was  used  on  vertical  concen¬ 
tration  data  and  the  results  used  to  estimate  depth  averages  and  surface-to- 
bottom  differences.  The  surface-to-bottom  difference  was  corrected  by  the 
inverse  of  the  fraction  of  the  water  column  that  was  sampled.  This  method  was 
checked  using  numerical  results  and  found  to  give  consistent  results.  The 
depth-averaged  diffusivity  for  homogeneous  flow  was  estimated  as 

K  =  0.067HU*  (26) 

z 

66.  CEC  is  a  measure  of  the  cations  held  on  the  surface  or  within  the 
crystaline  matrix  of  minerals.  It  is  affected  by  the  type  and  amount  of  clay 
minerals  present  as  well  as  some  organic  material.  Black's  method  of  analysis 
by  ammonium  saturation  was  used  because  it  minimizes  the  effects  of  organic 
material  on  CEC.  Analyses  were  done  by  the  Concrete  and  Materials  Analysis 
Group  of  the  Structures  Laboratory,  WES.  Typical  range  of  values  for  CEC 
found  in  the  literature  are 

_ Mineral  Fraction _  CEC,  meq/100  g 


Kaolinite 

1-15 

Illite 

10-40 

Montmorillonite 

50-150 

Chlorite 

10-40 

Vermiculite 

100-150 

Organic  fraction  of 

solids 

150-500 

Description 

of  Field  Methods 

and  Conditions 

67.  Field  sampling  for  sediment  characteristics  was  generally  carried 
out  in  conjunction  with  deployment  and  servicing  of  tide,  current,  and 
meteorologic  recording  instruments.  Sampling  was,  therefore,  neither  synoptic 
nor  time-series  in  nature.  Because  of  the  size  of  the  study  area,  usually 
more  than  1  day  was  required  to  complete  sampling.  Even  so,  not  every  station 
was  occupied  on  every  survey,  but  repeated  samplings  were  made  at  many. 
Sampling  methods  are  described  in  Coleman  et  al.  (1988).  Suspended  samples 
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were  generally  taken  with  a  pump  system  or  occasionally  with  a  Kemraer-type 
water  bottle.  On  one  occasion,  surface  grabs  were  taken  to  correlate  to 
satellite  imagery.  Fifteen  sampling  surveys  were  made  and  rated  somewhat  sub¬ 
jectively  according  to  the  transport  energies  imposed  by  diurnal  tide  range, 
wind,  freshwater  inflow,  and  TSM  (Table  7).  The  rating  covers  only  the  range 

of  observation  for  this  study.  For  instance,  high  inflow  in  the  table  is  con- 

3 

sidered  to  be  that  exceeding  6,400  m  /sec  (225,000  cfs)  while,  based  on  a 
longer  historical  record,  this  is  not  very  high.  River  inflow,  tidal  cur¬ 
rents,  and  waves  contribute  to  differing  degrees  to  the  TSM  levels  of  the  bay 
complex.  They  represent  sources  of  sediment  by  their  ability  to  transport  or 
erode.  For  comparison,  Figure  21  shows  a  relationship  between  freshwater 
flow,  Qj.  and  TSM  at  Morgan  City  while  Figure  22  shows  a  plot  of  freshwater 
flow  versus  TSM  for  river  mouth  stations  E  and  S.  The  latter  stations  were 
subject  to  more  tidal  effect  than  Morgan  City  and  therefore  can  be  expected  to 
have  greater  variability  in  the  discharge-concentration  relationship. 

Results  and  Discussion 

Tidal  influence  on  TSM 

68.  While  the  range  of  conditions  covered  by  these  observations  was 
limited  to  inflows  less  than  flood  level,  the  indications  are  that  sediment 
concentrations  in  the  bay  often  exceed  those  of  the  inflow.  This  is  common 
for  estuaries  in  general.  Plates  7-20  give  areal  distributions  of  TSM  for  the 
field  survey  sampling  periods  given  in  Table  7.  It  appears  that  tide  was  the 
most  dominant  single  factor  correlating  to  TSM,  although  inflow  discharge  cer¬ 
tainly  contributed.  Wind  had  surprisingly  little  effect  on  the  level  of  TSM 
in  the  bay.  The  most  severe  conditions  sampled  had  sustained  15-knot  winds 
accompanied  by  2-ft  waves.  Highest  concentrations  were  detected  at  times  of 
spring  tide  and  also  on  flood  phase  of  the  tide.  Figures  23  and  24  show  plots 
of  TSM  versus  diurnal  tide  range  at  stations  R  and  H,  respectively.  Flood  and 
ebb  phase  samplings  and  river  discharge  are  also  indicated.  Although  sampl¬ 
ings  were  not  taken  at  exactly  the  same  tidal  phase,  the  data  indicate  that 
flood  phase  concentrations  are  highest  and  that  an  upward  trend  of  concentra¬ 
tion  with  increasing  diurnal  range  exists.  The  correlation  of  TSM  to  river 
inflow  is  lower  at  these  stations  than  at  the  river  inflow.  These  stations 
are  very  much  tidal  at  low  discharges  and  data  are  scattered  accordingly. 
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Figure  22.  Freshwater  flow  versus  TSM  at  river  mouth 
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Figure  23.  TSM  versus  tidal  range  for  station  R 
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Figure  24.  TSM  versus  tidal  range  for  station  H 
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Distribution  of  settling  velocity 

69.  Table  8  summarizes  the  30  field  settling  tests  made  in  1980-1982. 
Figure  25  shows  the  areal  distribution  of  these  field  settling  velocity 
samples.  Figure  26  shows  a  plot  of  settling  velocity  versus  TSM  for  all  field 
settling  data.  The  correlation  of  settling  velocity  to  TSM  for  all  the  data 
is  very  poor;  but  at  individual  stations  the  correlation  might  be  better, 
although  not  enough  data  are  available  to  test  this  hypothesis.  In  general, 
settling  velocity  decreases  with  distance  away  from  the  river  mouths.  Fig¬ 
ure  27  shows  a  probability  plot  of  all  field  settling  data.  The  Weibull 
probability  distribution  displayed  the  data  in  a  more  straight-line  fashion 
than  did  normal  or  log-normr!  distributions.  Figure  28  shows  three  station 
settling  velocity  distributions  in  a  transect  out  from  Wax  Lake  Outlet  taken 
16  April  1982  that  demonstrate  this  trend.  Settling  ._iccity,  especially  in 
the  coarse  end  of  the  distribution,  and  concentration  both  fall  seaward  of  the 
river  mouth.  Settling  velocities  and  concentration  may  not  continue  to  fall 
into  the  offshore  area,  however.  Although  only  one  settling  velocity  test  was 
done  offshore  (station  L) ,  sediment  concentrations  and  stratification  were 
observed  to  be  high.  The  possibility  of  a  local  source  offshore  will  be 
further  discussed  in  the  next  section. 

Stratification  and  bed  shear 

70.  Part  II  presented  an  analysis  of  the  stratification  of  suspended 
material  in  terms  of  the  particle  Peclet  number  and  the  bottom  boundary  condi¬ 
tion  (erosion,  no  flux,  or  deposition).  Sufficient  data  were  collected  at 
some  stations  to  estimate  both  stratification  and  particle  Peclet  number; 
therefore,  some  inferences  as  to  the  sedimentation  conditions  at  the  bed  can 
be  drawn. 

71.  Figure  29  shows  a  diagram  for  the  classification  of  sediment  sus¬ 
pensions  for  these  stations.  Data  are  presented  in  Table  9.  Samplings  were 
taken  7-9  July  1982  at  midebb  phase  of  the  tide  with  2.0-  to  2.2-ft  diurnal 
range,  except  stations  U  and  V,  which  were  taken  near  high  water.  Freshwater 
inflow  was  225,000  cfs.  This  analysis  indicated  that  erosion  conditions  were 
occurring  at  all  stations.  Data  for  station  F  taken  11  May  1982  are  plotted 

for  comparison  and  indicate  a  flux  equilibrium  bottom  condition.  At  station  F 

2 

a  shear  stress  of  0.17  N/m  was  apparently  above  t  ,  the  c.  itical  shear 

2 

stress  for  deposition,  and  a  shear  stress  of  0.08  N/m  was  initiating  erosion 
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Figure  26.  TSM  versus  settling  velocity, 
all  field  tests 
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Figure  27.  Probability  distribution  of  field 
settling  data 


Figure  28.  Settling  velocity  distributions  for  three  stations 

across  the  hay 


at  station  U.  Thus  sediments,  at  least  more  mobile  fractions,  are  frequently 
in  transport  in  the  bay  system.  There  are  only  slight  differences  between  U 
and  U(sed)  in  these  data  so  that  the  assumption  that  sediment  is  moving  at 
the  depth-averaged  velocity  is  a  reasonable  one.  It  Is  apparent  that  at  least 
three  vertical  points  are  required  to  fit  vertical  profiles,  even  in  shallow 
water . 

72.  Table  10  gives  additional  data  on  shear  velocity  and  sediment 
stratification.  Offshore  stations  (L,  K,  and  tide  gage  8)  show  a  generally 
consistent  higher  sediment  stratification  than  inshore  stations.  At  the  times 
these  data  were  taken,  vertical  salinity  stratification  was  small,  less  than 
7  ppt  top  to  bottom  differences  at  station  L  and  tide  gage  8  in  25  ft  of  water 
and  even  less  at  station  K.  As  noted  in  the  last  section,  only  one  settling 
velocity  test  was  run  on  offshore  stations  that  indicated  a  high  settling 
velocity  with  mean  Wg  of  0.357  mm/sec  .  If  that  data  point  is  typical,  off¬ 
shore  stations  would  plot  on  the  erosional  side  of  the  suspension  classifica¬ 
tion  diagram,  indicating  an  offshore  sediment  source.  Offshore  sediment  beds 
are  commonly  (or  easily)  eroded,  providing  a  source  of  sediments  that  may  be 
then  transported  into  the  bays  for  deposition  or  transported  downcoast  and 
away  from  the  site. 
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Notes: 

Vertical  station  data,  collected  7- *9 Julu  198)2. 
All  mid  ebb  tide  samples  include  2.0'-  2.2'  ran^e 
except  stations  d  and  [/  (taken  near  raster'), 

t^restxrater  PJow  =  225^000  a/'s 


Figure  29.  Suspension  classification  diagram  for  a  high— transport  survey 


Bed  sediments 


73.  BWD  analyses  were  performed  on  vertically  sectioned  core  samples. 
Plates  21-28  give  the  vertical  profiles  of  BWD  at  the  stations  indicated  m 
Table  2.  CEC  data  are  given  in  Table  11.  Figure  30  is  a  plot  of  CEC  versus 
BWD  for  these  samples,  showing  a  tread  for  lower  CEC  and  higher  BWD  near  the 
river  mouths,  where  sediments  are  also  coarser.  High  CEC  and  low  BWD  were 
found  further  out  into  the  bay,  consistent  with  earlier  findings  that  the  fine 
fraction  is  most  cohesive  and  sediments  are  generally  fractionated  with 
coarser  material  near  the  river  mouths. 


Figure  30.  BWD  versus  CEC  for  Atchafalaya  sediments 


PART  V:  SUMMARY  AND  CONCLUSIONS 


74.  The  settling  and  consolidation  behavior  of  Atchafalaya  Bay  system 
bed  and  suspended  sediments  were  observed  under  various  conditions  in  labora¬ 
tory  and  field  experiments.  Four  ranges  of  concentration  affecting  settling 
were  identified.  The  effects  of  settling  height,  salinity,  and  turbulence  on 
settling  were  also  tested.  Settling  with  turbulence  and  hindered  settling 
consolidation  were  also  studied  using  numerical  models.  Measurements  in  bed 
sediments  were  made  for  bulk  wet  density,  mineral  size  distribution,  and 
cation  exchange  capacity.  Total  suspended  material  (TSM)  concentrations  were 
measured  in  suspended  samples. 

75.  The  following  conclusions  of  these  studies  as  they  apply  to 
Atchafalaya  Bay  were  drawn: 

a.  Bed  material  is  fractionated  in  the  bay  by  differential  set¬ 
tling  with  coarser  material  deposited  near  the  river  mouth  and 
finer,  more  cohesive  material  deposited  out  in  the  bay. 

b.  Settling  velocities  of  suspended  material,  in  general,  showed 
higher  rates  near  river  mouths  but  were  also  large  off  Eugene 
Island.  Field  settling  velocities  were  not  correlated  to  sus¬ 
pended  concentration  although  laboratory  tests  showed  that,  for 
a  given  sediment,  such  a  correlation  exists. 

c.  Over  the  conditions  sampled  in  the  field,  erosion  and  deposi¬ 
tion  are  cyclic  with  both  diurnal  and  spring-neap  tidal  cycles. 
Flood-tide  phase  sample  TSM  concentrations  are  generally  higher 
than  ebb  samples,  but  near  river  mouths  and  at  higher  river 
flows,  the  reverse  is  true.  Spring  tide  TSM  concentrations  are 
higher  than  neap.  Tidal  effects  were  found  to  be  stronger  than 
wind  or  even  river  discharge  over  the  range  of  conditions 
observed . 

d.  Mean  settling  velocities  of  suspended  material  were  found  to 
vary  from  0.005  to  0.3  mm/sec,  while  median  values  range  from 
0.001  to  0.1  mm/sec. 

e.  At  low  concentrations,  settling  is  independent  of  concentration 
(free  settling)  but  at  higher  concentrations  particles  inter¬ 
act,  increasing  aggregation  and  settling  (enhanced  settling). 

At  higher  concentrations,  particle  interactions  first  slow  set¬ 
tling  (hindered  settling)  and  eventually  make  particles  im¬ 
mobile  in  che  suspension  and  they  settle  as  a  mass  (hindered 
settling  consolidation).  Enhanced  settling  was  found  to  begin 
at  suspended  concentrations  about  0.1  to  0.2  g/£  and  continue 
through  2  or  3  g/£.  Hindered  settling  occurs  above  this  value 
and  was  found  to  follow  previously  established  empirical  laws 
(Equation  8).  Hindered  settling  consolidation  was  not  observed 
until  concentrations  of  the  order  of  100  g/ l  were  reached. 
Concentration  had  a  greater  effect  on  settling  velocity 
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(about  a  factor  of  10)  than  did  settling  height  (a  factor  of 
5),  salinity  (a  factor  of  2  to  A) ,  and  turbulence. 

2 

_f.  Turbulence  resulting  in  eddy  dif f usivities  of  285  mm  /sec  had 
no  marked  effect  on  the  equilibrium  settling  velocities  of  sus¬ 
pended  and  resuspended  bed  material.  Laboratory  and  field  set¬ 
tling  tests  are  therefore  felt  to  be  closely  related  to  in  situ 
equilibrium  settling,  although  field  settling  tests  tended  to 
give  results  about  a  factor  of  2  greater  than  laboratory  tests 
and  probably  are  the  most  realistic. 

£.  Field  estimates  of  the  particle  P  number  (W  H/K  )  range  from 
0.6  to  0.007.  Stratification  in  suspension  for  expected  inflow 
concentration  of  <1  g/i  is  expected  to  result  in  near-bottom 
concentrations  of  <1.5  g /£.  Since  the  onset  of  hindered  set¬ 
tling  consolidation  is  about  100  g/i.,  deposition  is  expected  to 
come  from  particle  settling  out  of  suspension. 

h.  For  the  observed  range  of  particle  and  suspension  stratif¬ 

ication,  sediment  transport  velocity  is  very  close  to  the  ver¬ 
tically  averaged  current  velocity;  and  the  free-settling  flux 
to  the  bed  approaches  an  exponential  function  of  depth  and  set¬ 
tling  velocity.  Therefore,  a  vertically  averaged  numerical 
modeling  approach  seems  valid. 

i.  Based  on  laboratory  tests,  initial  densities  of  newly^deposited 
sediment  are  expected  to  be  relatively  high  (0.7  g/cin  )  near 
river  mouths.  Out  in  the  bay  where  sediments  are  finer  and 
more  cohesive,  newly  deposited^sediment  densities  are  expected 
to  be  on  the  order  to  0.4  g/cm  . 

j_.  Hindered  settling  consolidation  of  thin  layers  of  sediment  with 
relatively  high  initial  concentrations  was  found  to  be  rapid. 
Numerical  results  showed  that  concentrations  can  double  in  only 
a  few  minutes. 
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Table  1 

Atchafalaya  River  Delta  Reports 


_ Title _ 

Plan  of  Investigation 

Field  Data  Program 


Extrapolation  of  Delta 
Growth 

Generic  Analysis  of 
Delta  Growth 


Quasi-Two-Dimensional 

Modeling 


laterim  Summary  of  Results 

Analytical  Analysis  of 
the  Development  of  the 
Atchafalaya  River  Delta 

Numerical  Modeling  of 
Hurricane-Induced  Storm 
Surge 

Wind  Climatology 
Wave  Hindcasts 
Two-Dimensional  Modeling 


2-D  Modeling  of 
Alternative  Plans  and 
Impacts  on  the 
Atchafalaya  Bay  and 
Terrebonne  Marshes 


_ Contents _ 

Investigation  Methods  and 
Approach 

Field  Data  Collection 
Methods  and  Presentation 
of  Data  (4  Sections) 

Analytical  Extrapolation 
of  Historical  Behavior 

Comparison  with  Similar 
Deltas  to  Identify  Stage 
of  Development  and  Predict 
Future  Trends 

Quasi-Two-Dimensional 
Hydrodynamic  and 
Sedimentation  Riverflow 
Modeling 

Summary  and  Analysis  of 
Reports  3,  4,  and  5 

Analytical  Treatment  of  a 
Simple  Jet  Discharging 
into  a  Quiet  Bay 

Two-Dimensional  Modeling 
Hurricane-Induced  Storm 
Surge 

Predictions  of  Wind 
Conditions  Over  the  Bay 

Modeling  of  Locally 
Generated  Wind  Waves 

Two-Dimensional  Finite 
Element  Modeling  of 
Hydrodynamics,  Salinity, 
and  Sedimentation 

Employs  the  tools  described 
in  Report  11,  and  shows  the 
effects  of  plans 
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Summary  Report 


Summary  and  Analysis  of 
the  Study 


Table  2 

Field  and  Laboratory  Sediment  Tests 
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X 
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*  BWD  *  Bulk  wet  density. 

**  CEC  =  Cation  exchange  capacity 
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Table  4 

Turbulence  Settling  Tests 


Screen 

Effective 

Initial 

Oscillation 

Settling 

Concentration 

Salinity 

Rate 

Velocity, 

mm/sec 

Test 

mg/ 1 

.  PPt 

Sta 

Source 

rpm 

Median 

Mean 

OA 

310 

0.2 

C,M 

Suspended 

90 

0.0085 

0.0196 

250 

0.0018 

0.0086 

440 

0.0016 

0.0045 

C3 

254 

4.3 

C,M 

Suspended 

90 

0.019 

0.0367 

250 

0.020 

0.0379 

440 

0.0049 

0.0152 

OC 

298 

3.8 

C,M,N 

Suspended 

90 

0.20 

0.0418 

250 

0.0078 

0.0204 

440 

0. 01)32 

0.0102 

OD 

177 

2.3 

C,M,N 

Suspended 

90 

0.0138 

0.0401 

250 

0.0073 

0.0315 

440 

0.0061 

0.0233 

OE 

198 

4.2 

M 

Bed 

90 

0.0369 

0.1652 

250 

0.0249 

0.1108 

440 

0.0170 

0.0847 

OF 

183 

0.2 

R 

Suspended 

0* 

0.0151 

0.0437 

90 

0.0080 

0.0389 

250 

0.0078 

0.0353 

440 

0.0047 

0.0224 

*  Performed  In  the  field  with  H  -  0.85  m 


Table  5 

Initial  Densities  of  Beds  Settled  from  Suspension 


Initial 

Concentration 

{?/* 

Settling  Height  Time 
cm  hr 

Bed  Thickness 

cm 

Bulk  Wet 
Density 

g  /  cm^ 

Concentration 

3 

g/cm  dry  wt . 

Station  M 

Composite 

78.8 

180 

2 

14 

1.52 

0.82 

39.5 

180 

4 

9 

1.45 

0.71 

78.8 

30.2 

0.2 

2.3 

1.45 

0.71 

39.5 

30.2 

0.1 

1.2 

1.42 

0.67 

8.2 

180 

0.5 

1.6 

1.46 

0.73 

41.1 

180 

0.5 

10.1 

1.41 

0.65 

Station  M 

Fine  Fraction 

3.8 

180 

0.5 

1.8 

1.13 

0.20 

3.4 

180 

0.5 

0.9 

1.31 

0.49 

27.8 

180 

0.5 

9.2 

1.30 

0.48 

Table  6 

Hindered  Settling  Consolidation  Tests 

Average 

Initial  Values 

Concentration 

BWD 

c 

- 

o 

H 

Duration 

C 

Station 

g/cm 

ill 

cm 

hr 

g  / 

T7 

X 

1.0980 

159 

160.0 

45 

o  ^  -» 

Z  /  / 

F 

1 .0887 

144 

100.0 

45 

533 

F 

1.0887 

144 

30.5 

50 

558 

H 

1.1100 

178 

160.0 

26 

380 

M 

1.1213 

202 

100.0 

28 

803 

M 

1.0864 

141 

100.0 

28 

742 

C 

1.0852 

139 

100.0 

54 

556 

C 

1.0852 

139 

30.5 

45 

695 

Table  8 


Field  Settling  Test  Conditions  and  Results  for  1980-1982 


Sta 

mg/ 1 

w 

s 

Cl 

C2 

C  3 

C  «  1 

SPD 

Dir 

Wind 

Sea 

A 

108 

0.1 

3 

17 

85 

0.2 

1.0 

090 

St  12 

— 

B 

13 

1.0 

3 

3 

7 

15.5 

0.0 

344 

00 

00 

C 

154 

0.0 

0 

20 

134 

0.0 

0.2 

323 

00 

0.2 

C 

126 

0.1 

5 

18 

103 

0.2 

0.4 

022 

SE  15 

2.0 

D 

39 

3.7 

1 1 

23 

5 

0.3 

0.5 

330 

NE  3 

0.3 

E 

33 

1.6 

3 

18 

12 

0.4 

1.0 

015 

N  2 

0.0 

F 

145 

3.5 

61 

33 

61 

0.2 

0.7 

040 

SE  15 

2.0 

F 

211 

2.3 

27 

138 

46 

0.0 

-- 

-- 

SE  10 

1.2 

H 

754 

3.0 

158 

453 

143 

0.3 

0.8 

042 

E  2 

00 

1 

198 

4.6 

59 

133 

6 

31.0 

0.2 

243 

00 

00 

K 

21 

1.6 

4 

8 

9 

31.0 

0.3 

170 

NE  4 

00 

K 

40 

0.7 

4 

14 

22 

1.3 

0.  7 

252 

SW  2 

0.8 

L 

198 

ir.0 

99 

91 

8 

19.3 

0.2 

314 

NW  6 

0.3 

M 

472 

1.5 

61 

222 

189 

0.7 

0.3 

098 

E  3 

0.3 

N 

154 

0.1 

15 

28 

111 

0.0 

0.4 

064 

SW  3 

0.8 

N 

93 

0.6 

7 

26 

60 

12.2 

— 

— 

E  15 

2.0 

C 

508 

4.9 

173 

310 

25 

0.0 

0.9 

052 

N  5 

0.8 

0 

159 

0.5 

14 

50 

95 

0.1 

-- 

SE  10 

1.0 

0 

100 

0.2 

8 

21 

71 

0.3 

0.2 

012 

NW  4 

0.7 

R 

392 

3.0 

90 

239 

63 

0.0 

1 .3 

330 

N  2 

00 

R 

560 

4.2 

168 

364 

28 

0.3 

0.7 

314 

SE  2 

0.2 

k 

500 

8.4 

235 

240 

25 

0.0 

— 

-- 

SE  10 

i  .  2 

R 

180 

i .  4 

29 

75 

76 

0.1 

— 

W  10 

1.0 

S 

12  3 

1 .  7 

18 

59 

46 

0.3 

0.5 

008 

00 

00 

s 

204 

0.4 

24 

5 1 

129 

0.0 

-- 

— 

SE  10 

— 

l1 

49 

1.9 

8 

24 

17 

9.4 

0.3 

014 

NW  3 

0.7 

f r 

1 1 : 

0.3 

20 

22 

69 

0.0 

0.4 

190 

W  1 

0.2 

V 

6  6 

0.4 

5 

23 

38 

0.0 

0.6 

1  18 

N  1 

0.5 

\ ' 

69 

0.3 

1  7 

1  1 

4  1 

5.9 

0 . 4 

303 

NW  3 

0.7 

X 

19 

0.2 

5 

3 

1  ! 

15.0 

0 . 6 

290 

jE  10 

1.0 

Note:  Cl  =  TSM  with  settling  velocity  greater  than  10  mnw'sec  (equivalent 

Stokes  diameter  >  12  um)  ;  , 

C2  =  TSM  with  settling  velocity  between  10  and  10  ~  mm/see  (equiva¬ 
lent  Stokes  diameter  between  12  and  3,5  urn’) ; 

_ 

C3  =  TSM  witn  .settling  velocity  ’ess  than  10  mai/sec  (equivalent 
Stokes  diameter  <  3.5  urn); 

Sal  -  salinity,  ppt ; 

SP1J  =  current  speed,  fps; 

Dir  =  current  direction,  deg; 

Wind  direction  and  snnec  iti  knots; 

Sea  -  wind  waves ,  ft. 


Table  9 

Suspension  Classification  Data  for  Stations  Diagrammed  in  Figure  29 


Station 

Date 

Depth 

m 

U 

ra/sec 

C 

mg/  £ 

U*(fm) 
m/ sec 

Stratification 

P 

e 

F 

7/07/82 

5.34 

0.448 

322 

0.021 

1.53 

0.072 

R 

7/07/82 

1  nrt 

1  *  O  U 

0.561 

371 

0.032 

A  A  » 

>J  •  >  M 

yj  •  C  u 

0 

7/07/82 

1.98 

0.282 

433 

0.016 

2.80 

0.171 

H 

7/07/82 

2.93 

0.254 

379 

0.013 

1.07 

0. 101 

V 

7/09/82 

2.29 

0.206 

55 

0.011 

0.90 

0.024 

U 

7/09/82 

2.90 

0.163 

112 

0.009 

0.17 

0.025 

U 

1/25/83 

2.29 

0.270 

202 

0.015 

0.51 

0.220 

U 

2/23/83 

2.59 

0.035 

536 

0.002 

0.64 

0.007 

R 

2/24/83 

1.83 

0.426 

183 

0.024 

0.48 

0.010 

0 

2/24/83 

2.44 

0.090 

162 

0.005 

1.25 

0.242 

H 

3/29/83 

2.44 

0.304 

185 

0.016 

1.00 

0.119 

H 

4/26/83 

2.44 

0.419 

277 

0.023 

0.68 

0.086 

0 

4/26/83 

1.83 

0.093 

93 

0.005 

0.27 

0.012 

R 

5/24/83 

1.98 

0.874 

255 

0.049 

0.A3 

0.043 

H 

5/24/83 

2.44 

0.487 

314 

0.026 

0.16 

0.073 

C 

5/24/83 

1.22 

0.206 

180 

0.012 

0.71 

0.016 

H 

6/22/83 

3.66 

0.516 

167 

0.026 

0.76 

0.100 

R 

6/22/83 

2.44 

0.615 

141 

0.033 

0.42 

0.074 

r 

6/22/83 

1.83 

0.244 

76 

0.014 

0.27 

0.084 

r 


Table  10 

Shear  Velocity  and  Suspension  Stratification  Data 


Station 

Date 

Depth 

m 

U 

m  /  sec 

C 

g/ 1 

6C/C 

A 

05/11/82 

2.0 

0.336 

0.019 

0.  Iu3 

0.0 

C 

01/24/83 

2.3 

0.248 

0.014 

0.210 

0.701 

C 

05/24/83 

1.2 

0.206 

0.022 

0.180 

0.712 

C 

06/22/83 

1.8 

0.244 

0.014 

0.076 

0.265 

F 

05/11/82 

5.3 

0.266 

0.018 

0.109 

0.972 

F 

06/08/82 

4.9 

0.223 

0.013 

0.326 

1.064 

F 

07/07/82 

5.3 

0.448 

0.021 

0.322 

1.530 

H 

08/05/80 

2.4 

0.417 

0.023 

0.033 

0.28 

H 

05/11/82 

2.7 

0.473 

0.025 

0.077 

0.32 

H 

06/08/82 

2.9 

0.224 

0.012 

0.221 

0.13 

H 

07/07/82 

2.9 

0.254 

0.013 

0.379 

1.070 

H 

02/24/83 

2.4 

0.270 

0.015 

0.194 

0.229 

H 

03/29/83 

2.4 

0.304 

0.016 

0.185 

1.004 

H 

04/26/83 

2.4 

0.419 

0.034 

0.277 

0.677 

H 

05/24/83 

2.4 

0.487 

0.024 

0.314 

0.15? 

H 

06/22/83 

3.7 

0.516 

0.043 

0.167 

0.764 

I 

05/11/82 

1.7 

0.183 

0.011 

0.072 

0.18 

T 

06/09/82 

1.7 

0.137 

0.008 

0.312 

0.06 

K 

06/09/82 

3.4 

0.112 

0.005 

0.130 

1.38 

K 

07/07/82 

3.1 

0.143 

0.007 

0.061 

1.06 

K 

04/26/83 

2.0 

0.286 

0.016 

0.104 

n .  i  5 

L 

08/05/80 

6.4 

0.447 

0.021 

0.042 

2.51 

L 

08/19/80 

7.0 

0.16G 

0.008 

0.028 

1  .71 

L 

06/09/82 

6.9 

0.163 

0.009 

0.015 

2.59 

I. 

07/07/82 

7.0 

0.203 

0.009 

0.010 

3.54 

M 

06/08/82 

1 . 7 

0.107 

0.006 

0.151 

0.36 

M 

«.  '08/8  2 

1.6 

0.092 

0.005 

0.  159 

NA 

N 

06/09/82 

2.0 

0.102 

0.005 

0.244 

0.397 

0 

08/05/80 

1.5 

0.488 

0.029 

0.095 

0.66 

0 

08/19/80 

1  .5 

0.279 

0.013 

0.022 

0.36 

n 

06/10/82 

1.8 

0.163 

0.0C9 

0.34  5 

1.82 

(Continued) 


Table  10  (Concluded) 


Station 

Date 

Depth 

m 

U 

m/ sec 

U' 

m, 

EfliS 

C 

g/£ 

6C/C 

0 

07/07/82 

2.0 

0.282 

VJ 

.016 

0.433 

2.801 

C 

01/24/83 

2.8 

0.122 

0 

.006 

0.085 

0.710 

0 

02/24/83 

2.4 

0.090 

0 

.005 

0. 162 

1.252 

0 

04/26/83 

1.8 

0.093 

0 

.005 

0.093 

0.268 

0 

05/24/83 

2.4 

0.284 

0 

.015 

0.102 

0.936 

R 

08/05/80 

1.5 

0.549 

0 

.032 

0.171 

0.30 

R 

08/19/80 

1.6 

0.336 

0 

.020 

0.062 

0.98 

R 

06/09/82 

1.8 

0.254 

0 

.015 

0.903 

0.71 

R 

07/07/82 

1.9 

0.561 

0 

.032 

0.371 

0-943 

R 

02/24/83 

1.8 

0.426 

0 

.024 

0.183 

0.482 

R 

05/24/83 

2.0 

0.874 

0 

.049 

0.255 

0.432 

R 

06/22/83 

2.4 

0.615 

0 

.033 

0.141 

0.416 

U 

08/06  '80 

2.9 

0.203 

0 

.011 

0.042 

0.49 

u 

06/09/82 

3.1 

0.112 

0 

.005 

0.059 

0.27 

u 

07/09/82 

2.9 

0.163 

0 

.009 

0.112 

0.171 

rj 

01/25/83 

2.3 

0.270 

0 

.015 

0.202 

0.505 

u 

02/23/83 

2.6 

0.035 

0 

.002 

0.536 

0.636 

V 

08/06/80 

2.9 

0.163 

0 

.008 

0.045 

0.12 

V 

07/09/82 

2.3 

0.206 

0 

.01  1 

0.055 

0.902 

V 

01/25/83 

2.13 

0.094 

0 

.005 

0.225 

1  .455 

V 

04/27/83 

2.3 

0.132 

0 

.007 

0.147 

0.002 

X 

08/06/80 

2.9 

0.427 

0 

.022 

0.053 

0.20 

Tide  Gage  8 

08/08/80 

8.8 

0.366 

0 

.016 

0.055 

2.39 

Y 

08/06/80 

1.8 

0.671 

0 

/■* 

0.033 

0.32 

Table  11 

Test  Results  for  Cation  Exchange  Capacity 


Sample 

Station 

Depth  in  Sediment 
ft 

CEC 

meq/100  g 

F 

Composite 

27. A 

F 

0.1 

35.7 

F 

0.6 

21.9 

R 

0.1 

14.9 

R 

0.6 

15.0 

G 

0.1 

17.2 

G 

0.6 

10.0 

X 

0.75 

5.8 

K 

Composite 

35.4 

K 

0.1 

38.7 

K 

0.6 

24.7 

M 

Composite 

11.2 

M 

0.1 

10.5 

M 

0.6 

8.0 

H 

Composite 

31.6 

C 

Composite 

42.4 

C 

0.1 

20.7 

C 

0,6 

32.2 

Commercial 

(Kaolinite) 

13.1 

13.1 


Station  F 


Comparison  Test  T 
Station  F 


Comparison  Test  S 
Station  F 


Comparison  Test  T 
Station  C 


PLATE  6 


Comparison  Test  S 
Station  C 


PLATE  8 


PLATE  10 


27- 29  Ocf  /?«/ 

-  refer  fo  Ft  a  arc  £  -far 


TSM,  27-29  Oct  1981 


PLATE  12 


station  designations 


TSM,  11  May  1982 


o 


-0.4  - 


$  -0.& 


Dry  density  ,  !t>  dry  neuyht/fy.3 

20  40  40  30  /< 

_ 1 _ 1 _ 1 _ ! _ 

Bulk  epee/ fie  cynav/ty  2 


soft  sandy  day 

[bronrr  and* yrtny') 

ft  sandy  dau 
( branny  ’ 


Soft  sandy  da 
( branny 


soft  day  nith  some  aha/I 
Cyrdy) 

medium  stiff  ctau  nith 
Some  shall  (ylray) 

stiff  ctau  wth  sha/is 
(a ratf) 


Q  '  °a\ 


Stiff  day 

l  yrayi) 


^  -ao 

< 

'V 

S 


STATION  A 


Dry  density  ,  lb  dry  neiyht/fyi 

20  40  40  BO  ZOO 


Bulk  specific  ynautty 


Soft  day,  htt/a  sand 

(J*ranrif 

Sti  ff  day  nith  sand 
(brand) 

stiff  day  wth  »and 


rpjtfn  m 

,  fray) 


stiff  da, 


iff  day 
[brand,  yray) 


STATION  C 


BULK  WET  DENSITY 
Stations  A  and  C 


PLATE  21 


Dry  density  ,  lb>  dry  we/yht^f+S 


S  -o.&\ 


§-aa 


Bulk  specific  gravity 


STATION  H 


c/ay,  ssffi*  shell 
(  brown ,  yro  y  ) 

Stiff  c/eu ,  some  shell 
( brown,  tfroy) 

\  oarj ,  some  shell 
/yray) 

c>4iy  »v//A  shell 

'(fay) 

clou  sand  with  shell 


Dry  density  ,  /A  w/e/yht^f^ 

O  20  40  <kO  SO  /OO 


Bulk  specific  cynovity 


soft-  dau  , 

C  brow/1,  yray) 

with  sand 


-C 

'is 


STATION  I 


BULK  WET  DENSITY 
Stations  H  and  I 


PLATE  23 


PLATE  2k 


Dry  density  ,  ib  dry  tuetyht^f+3 

20  40  60  1  ao  too 


Bulk  specific  eyrsvity  Z.O 


brotvn  0rau  soft  den. 


bromn  qrau  stiff  clou 

h*He< sathd  7 


bnsivn  arau  stiff  da. 


Sandy  day 


STATION  N 


Dry  density  ,  lb  dry  tueiyfit/f+b 

20  40  60  BO  /OO 

I _ I _ I _ I _ I 

Bulk  specific  yrtx  v/ty  2-0 

brorjn  sandy 

^  bromn  sandy 

/  bron/n  sandy  day 


STATION  0 


BULK  WET  DENSITY 
Stations  N  and  O 


PLATE  25 


Ory  density  ,  /b.  dry  ivet'yhl foa 

O  20  40  60  go  /CO 

1 _ ! _ I _ '  j 


1.0 

o\— 


0.2 


-0.4 


£  -06 


-  1.0 


2.0 


soft  cJau  *v///>  sand 
(broivA) 

sand  and  c/au 

(yray  ana  rbrotvn) 

Soft  claq  *v/M  sand 
(,  1r*y  ‘and  bro*jn  ) 


Sandy  day 

dayey  sand 

(<?ray  and  brown) 


clayey  sand 

(•fray  and  hrown  'j 


STATION  O 


Dry  density  ,  Id  dry  tvetyht^^a 

O  20  40  60  go  /OO 


1.0 

O  h 


-0.2 


-0.4 


20 

_J 


dark  tiroivn 
she/I  and  sand 

sand 

bnotvn  anau  sand 
*v/M  SkeV/ 

bra  tun  ycay  sand 


^  -06 


-  i-0 L  STATION  R 


PLATE  26 


BULK  WET  DENSITY 
Stations  O  and  R 


Dry  density  ,  lb  dry  meiyhi/f+i 


Bulk  spacing:  cynxvity 


ck>u  and  sard  with  aome 
'  shell  ( brand) 


clay  (brawn) 

c/om  with  soma  sand  and 
'  shell  (yray) 


yray  day  with  soma  shell 


STATION  R 


Dry  density  ,  Jb  dry  weiytt 

20  40  60  '  30  IOO 


Bulk  specific  eynavity 


clayey  aitt 

SiH  with  some  cJay 

Silt 

silt  witty  some  day 
c/ayey  aitt 

clayey  si  It 


STATION  S 


BULK  WET  DENSITY 
Stations  R  and  S 


PLATE  27 


PLATE  28 


BULK  WET  DENSITY 
Stations  U  and  *NEWa 


APPENDIX  A:  NOTATION 


Subscripts  and  indices: 
b  Bottom;  bed 

f  Fresh  water 

o  Indicates  a  depth  average  or  vertically  homogeneous 

x  Horizontal  coordinate 

y  Vertical  coordinate 

z  Depth  coordinate 

Parameters : 

a  Regression  constant 

BWD  Bulk  wet  density 

C  Concentration  of  suspended  sediment 

C  Depth-averaged  concentration 

C  Initial  concentration 

o 

C  Concentration  of  bed  sediments 

s 

CEC  Cation  exchange  capacity 

D  Deposition 

E  Erosion 

g  Force  due  to  gravity 

h  Deviation  from  the  mean  tide  level 

H  Thickness  of  mud  layer;  total  depth  or  height;  settling  height 

k  An  index,  such  as  the  number  of  sediment  fractions 

K  Empirical  coefficient 

Effective  vertical  eddy  diffusivity 
m  Mass  rate  of  erosion 

m  Mass  per  unit  surface  area  of  the  consolidating  layer 


Ai 


n 

P 


P 

e 

S 


t 

T 

TSM 


U 

u* 

u 


U*(fm) 


W 

s 


W 

s 

g 

W 

s 

X 

X 

z 

Z 


Z 

o 

6C 

P 


P 

w 


P 

o 


a 


a 

g 


Manning's  friction  coefficient 
Probability  of  particle  sticking  to  bed 
Particle  Peclet  number 
Freshwater  inflow 
Salinity 

Nondimens ional  time 
Time,  dimensional 
Total  suspended  material 
Horizontal  velocity 
Bed  shear  velocity 
Depth-averaged  velocity 

Shear  velocity  estimated  using  Manning's  equation 
Settling  velocity 
Hindered  settling  velocity 
Geometric  mean 

t  ll 

Settling  velocity  at  the  x  percentile 

Nondimensional  horizontal  coordinate 
Nondimensional  vertical  coordinate 
Vertical  dimension  measured  upward  from  bed 
Roughness  height 

Surf ace-to-bottom  difference  in  concentration 
Density 

Density  of  water 

Reference  density  (fresh  water) 

Dye  variance 


Geometric  standard  deviation 


Shear  stress 


Bed  shear  stress 

Critical  shear  stress  for  deposition 


